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ABSTRACT 

We compute rates of tidal disruption of stars by supermassive black holes in galactic nuclei, using downwardly 
revised black hole masses from the M, - a relation. In galaxies with steep nuclear density profiles, which dominate 
the overall event rate, the disruption frequency varies inversely with assumed black hole mass. We compute a total 
rate for non-dwarf galaxies of ~ 10~ 5 yr -1 Mpc" 3 , about a factor ten higher than in earlier studies. Disruption 
rates are predicted to be highest in nucleated dwarf galaxies, assuming that such galaxies contain black holes. 
Monitoring of a rich galaxy cluster for a few years could rule out the existence of intermediate mass black holes 
in dwarf galaxies. 

Keywords: galaxies: elliptical and lenticular — galaxies: structure — galaxies: nuclei — stellar dynamics 



1. INTRODUCTION 

Stars that pass sufficiently close to a supermassive black 
hole will be tidally disru pted (Hi lls 19751 iFrank & Rees 19761 
Lidskii & Ozernoi 19791 . Disruption of solar-type stars oc- 
curs at a distance r, « ^(M./OVio) 1 / 3 , with M. the black 
hole mass; for M, < 10 8 M Q , the tidal radius lies beyond the 
black hole's event horizon and disruption results in an ener- 
getic flare as the bound stellar debris falls back onto the black 
hole. Emission from the debris is expected to peak in the 
soft X-ray or UV domains, to have a maximum luminosity of 



to months JRees 19881 


Evans & Kochanek 1989 Ulmer 1999 


|Kim, Park & Lee 1999] 


. Detection of flares would constitute 



robust proof of the existence of supermassive black holes and 
could conceivably allow constraints to be placed on black hole 
masses and spins (Rees 1998 1. 

The ROSAT All-Sky Survey detected soft X-ray outbursts 
from a number of galaxies with no previous history of of Seyfert 
activity. Roughly half a dozen of these events had the proper- 
ties of a tidal disruption flare (Komossa 2002 and references 
therein), and follow-up optical spectroscopy of the candidate 
galaxies confirmed that at least two were subsequently inactive 
dGezari et al. 2003t . 

The mean event rate inferred from these outbursts is roughly 
consistent with theoretical predictions (Donley et al. 2002). 
Detailed calculations of the tidal disruption rate in samples of 
nearby galaxies have been published by Syer & Ulmer (1999, 
hereafter SU) and Magorrian & Tremaine (1999, hereafter 
MT). Both groups took black hole masses from the Magor- 
rian et al. (1998) demographic study, which found a mean ra- 
tio of black hole mass to bulge mass of ~ 0.006. Following 
the discovery of the M, - a relation ( Ferrarese & Merritt 2000 
Geb hardt et al. 20001 . the mean value of M,/Mb u ige was 
revised downward, to ~ 0.001 dMerritt & Ferrarese 20 01a 
Kormendy & Gebhardt 2001 1. The lower mean value of 
M,/Mbuige resolved two outstanding discrepancies: the factor 
~ 10 difference between black hole masses in quiescent and ac- 
tive galaxies with similar luminosities (Wandel 1999); and the 
higher apparent density of black holes in nearby galaxies com- 
pared with what is needed to explain the integrated light from 
quasars (Rich stone et al. 1998i . 



Here we examine the consequences of downwardly revised 
black hole masses for the rate of stellar tidal disruptions in 
galactic nuclei. Published scaling relations (e.g. Frank & Rees 
1976; Cohn & Kulsrud 1978) suggest that stellar consumption 
rates in galactic nuclei should scale as ~ MJj,4/3 < n < 9/4 
when the other properties (density, velocity dispersion) of the 
host galaxy are fixed. Hence one might naively expect the lower 
values of M. to imply lower rates of stellar disruption. Instead 
we find the opposite: in most galaxies, and in particular in the 
galaxies with steep nuclear density profiles that dominate the 
overall event rate, decreasing the assumed value of M m leads to 
higher rates of loss cone feeding. We estimate a total tidal dis- 
ruption rate among non-dwarf galaxies that is about an order of 
magnitude higher than in studies based on the Magorrian et al. 
black hole masses. 

This paper is laid out as follows. §2 describes the galaxy 
sample and §3 reviews the steady-state loss cone theory from 
which event rates are computed. The theory is applied in §4, 
with the counter-intuitive result that lower values of M. imply 
greater feeding rates in most galaxies. This result is analyzed 
in more detail in §5, where it is shown to be a generic property 
of steep power-law nuclei. We derive an accurate, analytical 
expression for the tidal disruption rate in singular isothermal 
sphere nuclei. In §6 we present the implications of our results 
for the overall rate of tidal flaring and show that the predicted 
rate would be so high in dwarf galaxies that the presence of 
black holes in these systems could be ruled out by just a few 
years' monitoring of a rich galaxy cluster like Virgo. §7 sums 
up. 

2. GALAXY SAMPLE 

Our basic sample is the set of 61 elliptical galaxies whose 
surface brightness profiles were studied by Faber et al. (1997). 
These authors fit the luminosity data to the parametric model 



(i) 



where rj, is the "break radius," 4 = /(r/,), and V is the logarith- 
mic slope of the surface brightness profile at small radii. (Note 
that we adopt the "theorist's convention" in which 7 refers to 
the logarithmic slope of the central space density profile.) For 
51 of these galaxies, Faber et al. (1997) give values for each 
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of the five parameters /i£,a, (3, Y and Yy; the latter is the vi- 
sual mass-to-light ratio assuming Hq = 80 km s" 1 Mpc" 1 and 
fib is the surface brightness at rb in visual magnitudes arcsec -2 . 
For these 5 1 galaxies we computed the mass density profile via 
Abel's equation: 



p(r) = T v j(r) = - 



T, 



dl dR 



(2) 



with j(r) the luminosity density. Below we follow the con- 
vention of referring to galaxies with T < 0.2 as "core" galax- 
ies and those with Y > 0.2 as "power-law" galaxies. 1 We 
note that even most core galaxies exhibit an approximately 
power-law dependence of space density on radius for small r, 
p ~ r~'< (Mer ritt & Fridman 19961 IGebhardt et al. 19961 . The 
weak power-law dependence p on r in the core galaxies is 
not well reproduced by deprojection of the fitting function Q, 
which is a possible source of systematic error in what follows. 
Our sample (Table 1) contains 28 power-law galaxies and 23 
core galaxies. 

The gravitational potential ip(r) = — $(»*) was computed via 



^(r) = ^*(r) + 



GM. 



(3a) 



4ttG 



p(r')r' dr+AixG / p(r'ydr' (3b) 

r dr 



GM(r) A „ f°°T v f°°dI(R) dR , , 



-4nG I — 

7T 



dR V/? 2 -r' 2 



(3c) 



a \ / - 7 - l 



K + r'y-^W^-^dr' 



(3d) 



The distribution function /, defined as the number density of 
stars in phase space, was computed via Eddington's formula: 



/(e) = 



1 



dp dip 



%/87r 2 m* de J dtp yje-ip 



(4) 



with m* the stellar mass and s = —E. We follow MT in assum- 
ing an isotropic velocity distribution. Galaxies with Y < 0.05 
were found to have / < when M, > 0; this is a consequence 
of the fact that an isotropic / can not reproduce a shallow den- 
sity profile around a point mass. The 10 galaxies with negative 
f's are included at the end of Table 1 and not discussed further 
here. 

We define the sphere of influence of the black hole to have 
radius rn, where 

MJr h ) = 2M. (5) 



A second way to estimate black hole masses is via the M m — a 
relation (Ferrarese & Merritt 2000; Gebhardt et al. 2000). We 
adopt the updated version from the review of Merritt & Fer- 
rarese (2001b): 



M, = 1.48 x 10 M.c 



200 km s- 



(7) 



the errors in the normalizing coefficient and exponent are 
±0.24 x 10 8 .Mq and ±0.48 respectively. Equation Q was de- 
termined from a fit to the small sample of galaxies in which 
the black hole's sphere of influence is clearly resolved. The pa- 
rameter cr c in equation (Q is the velocity dispersion measured 
in an aperture of si ze r e /8 centered on the nucleus, with r e the 
effective radius (Ferrarese & Merritt 2000). We computed a c 
from published measurements of the central velocity dispersion 
following the prescription in Ferrarese & Merritt (2000). The 
M. — a masses are listed in column 13 of Table 1. 

As can be seen in Table 1, and discussed in detail elsewhere 
(e.g. Ferrarese & Merritt 2000; Merritt & Ferrarese 2001b), the 
Magorrian et al. masses are systematically high compared with 
masses computed via the M. - a relation. It is this discrepancy 
which motivated the current study; both published studies of 
stellar disruption rates in galactic nuclei (SU, MT) were based 
on the Magorrian et al. masses. 

3. LOSS CONE THEORY 

Stars of mass m* and radius r* that come within a distance 



(j] ) 7* 

rn,. 



(8) 



of the black hole will be tidally disrupted; r\ w 0.844 for an n = 3 
polytrope. Following MT, we define the "consumption rate" 
as the rate at which stars come within r,, even if r, falls be- 
low the Schwarzschild radius 2GM,/c 2 ; the latter occurs when 
M, > lO 8 ^©. (The largest consumption rates occur in small 
dense galaxies for which r, > r s .) In a spherical galaxy, stellar 
orbits lie within the consumption loss cone if their energy e and 
angular momentum per unit mass J satisfy 



J 2 < 4(e) ^ 



2rf[iP{r t )-e}~2GM.r,. 



(9) 



We adopt the Cohn-Kulsrud (1978; hereafter CK) formalism 
for computing the flux of stars into the loss cone. Let T{e) 
be the number of stars per unit time and unit energy that are 
deflected into the loss cone via gravitational encounters with 
other stars. Define ((AR) 2 } to be the diffusion coefficient in 
R = J 2 /J 2 (e), with / c (e) the angular momentum of a circular 
orbit of energy e. Then (CK) 



andM + is the mass in stars within r. This definition is equivalent 
to r/, = GM,/a 2 when p(r) oc r~ 2 . We further define £/, = ip(rh). 

Of the 41 galaxies in our sample with non-negative f's, 18 
have black hole masses tabulated in Magorrian et al. (1998). 
These masses are given in column 1 1 of Table 1 . For the re- 
maining galaxies, we give in column 11 black holes masses 
computed from 

M. = 0.006M hu , ge (6) 

the mean relation between bulge mass and black hole mass 
found by Magorrian et al. (1998). 

'in fact we retain the Faber et al. (1997) classifications in Table 1, which are slightly different. 



f(e)de = 4ir 2 J 2 



dr ((AR) 2 ) 

— hm 

V r R-^o 2R 



f 



Info 



-de 



and the total consumption rate is given by 



7V = 



T(e)de. 



(10) 



(ID 



In equation (1101 . Rq(e) is the value of R at which / falls to zero 
due to removal of stars that scatter into the loss cone. Rq is not 
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equal to its geometrical value, Ri c = jf c /J c (e) 2 , because scatter- 
ing of stars onto loss cone orbits permits / to be nonzero even 
for J < Ji c . CK find that Ro(e) can be well approximated by 



R Q (e) = R lc (e)x 



f exp(- 
\ exp(- 



q), 

-0.186?- 



0.824^), 



q(£) > 1 
q(£) < 1 



with 



1 



Ri c (e) J v r r^o 2R R, c (e) 



(12) 



(13) 



P(e) is the period of a radial orbit with energy e and ~p is the 
orbit-averaged diffusion coefficient. The function q(e) can be 
interpreted as the ratio of the orbital period to the time scale for 
diffusional refilling of the loss cone; q ^S> 1 defines the "pin- 
hole" or "full loss cone" regime in which encounters replenish 
loss cone orbits much more rapidly than they are depleted. MT 
give expressions for the local angular-momentum diffusion co- 
efficient: 

((AR) 2 ) 32n 2 r 2 G 2 m 2 JnA 



/o 



lim 

R^o 2R 



f(e')de', 



3^ (3/ 1/2 -/ 3/2 + 2/o), 

(14a) 
(14b) 



[2(^(r)- £ ')] 3/2 /(£')^' (14c) 

from which the orbit-averaged quanties 71(e) and q(e) can be 
computed. 

In the Fokker-Planck approximation under which equation 
dipt was derived, the flux of stars into the loss cone at each en- 
ergy is determined by gradients in / with respect to R at R « 
Ric- CK, who modelled globular clusters, derived expressions 
for these gradients by assuming that the distribution of stars 
near the loss cone had evolved to a steady state in which the 
encounter-driven supply of stars into the loss cone was balanced 
by consumption. Relaxation times in galactic nuclei are usu- 
ally in excess of a Hubble time, particularly in the core galaxies 
( Fab er et al. 19 97 1, and it is not clear that the distribution of 
stars near the loss cone will have had time to reach a steady 
state in all of our galaxies (Milosavljevic & Merritt 2003 1. We 
will return to this question in a subsequent paper; for now, we 
follow MT in assuming that the CK loss cone boundary solution 
applies to galactic nuclei. 

4. DEPENDENCE OF THE CONSUMPTION RATE ON M. IN THE 
GALAXY SAMPLE 

Figures 1 and 2 show how the energy dependence of various 
quantities changes with the assumed value of M. in two galax- 
ies: NGC 4551, a power-law galaxy (r = 0.8); and NGC 4168, 
a core galaxy (T = 0. 14). The value of Yy , and hence the mass 
density of the stars, was fixed as M m was varied. We adopted 
m+ = A4q throughout. In the power-law galaxy, as M, is de- 
creased, the flux of stars into the loss cone increases at e > £/, 
and decreases at e < £/,; since most of the flux comes from near 
the black hole, £>£/,, the total consumption rate increases with 
decreasing M, . In the core galaxy, the dependence of T{e) on 
M. is more complex: T(e ~ £/,) first increases, then decreases, 
with decreasing M,. The consequences of these trends can be 
seen in Figure 3, which plots integrated consumption rates N 





1 




log (£/(100kms 1 ) 2 ) 



log (£/(100kms 1 ) 2 ) 



FIG. 1. — Dependence on M, of various quantities associated with stellar 
consumption in the power-law galaxy NGC 4551. Stars are assumed to have 
solar mass and radius. Left column: r apo (apocenter radius of radial orbit); 
/ (phase space number density); q (quantity that distinguishes between diffu- 
sion and full-loss-cone regimes). Right column: P (period of radial orbit); Ri c 
(geometric size of loss cone in terms of R = J 2 /J^); T (flux into loss cone). 
As M, is reduced, more and more of the galaxy falls within the full-loss-cone 
regime (q ^> 1) and the total flux of stars into the loss cone rises. 
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FIG. 2. — Like FigureQ but for the core galaxy NGC 4168. By comparison 
with NGC 4551, less of the galaxy lies in the full-loss-cone regime and the 
total consumption rate is lower. 
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as a function of M, for every galaxy in our sample. The power- 
law galaxies exhibit monotonic or nearly montonic trends of 
increasing N with decreasing M.; the dependence is roughly 
N oc M~ l . In the case of the core galaxies, N generally increases 
with M, up to a maximum value, then decreases as M, is in- 
creased further. This behavior is explained in the next section. 
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Fig. 3. — Dependence of consumption rate on assumed black hole mass for 
the galaxies in Table 1 . 

Our primary concern is how consumption rates would change 
if the black hole masses used by earlier authors were replaced 
with the presumably more accurate masses derived from the 
M, - a relation. Figure 4 makes this comparison. In almost 
every galaxy in our sample, the inferred N is greater when the 
M, - a black hole mass is used. The changes are greatest in 
the power-law galaxies, since N ~ M~ l in these galaxies. In 
the core galaxies, although M, is sometimes increased by as 
much as 10 2 , the changes in N are usually modest because of 
the nearly flat dependence of N on M, in these galaxies (Figure 
3). 



Figure [5] shows the dependence of N on galaxy luminosity 
and M. in our sample; black hole masses were taken from the 
M. - cr relation. The dependence of flaring rate on M. is fairly 
tight, with a mean slope of N ^ M~ 0,8 . 
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FIG . 5 . — Consumption rate as a function of galaxy luminosity (a) and black 
hole mass (b). The dashed line in (b) is the relation defined by the singular 
isothermal sphere, eq. (29); it is a good fit to the galaxies plotted with stars, 
which have steep central density profiles, p ~ r~ 2 . 
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log 10 N (yr 1 ) (M. - a) 

FIG. 4. — Comparison of consumption rates computed using the two values 
of M. in Table 1 . Abscissa: M. computed from the M, - a relation, equation 
[7J. Ordinate: M. from Magorrian et al. (1998). 



Figure[6]shows how and three critical radii associated with 
the black hole vary with M. in NGC 4551 and NGC 4168. The 
tidal radius r, and the radius of influence r/, were defined above. 
The third radius, r crit , is defined as: 



ip(r cr it) = £c 



q(Scrit) = 1 • 



(15) 



r cr i, is roughly the radius of transition between the "diffusive" 
(q < 1) and "full-loss-cone" (q > 1) regimes, and most of the 
flux into the loss cone comes from radii r < r cr it- We note that 
r crit < r/, over the relevant range in M, for both galaxies, and 
that both r, and r crl , are less than n,. These same inequali- 
ties were found to hold for most of the galaxies in our sam- 
ple, which motivated the simplified treatment in the following 
section. 

5. DEPENDENCE OF CONSUMPTION RATE ON M. IN POWER-LAW 

NUCLEI 

We noted above the curious behavior of as M. is varied in 
a galaxy with otherwise fixed properties: generally increases 
as M, is reduced. Here we show how the dependence of on 
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M. can be understood. We derive the exact consumption rate in 
a p oc r~ 2 galaxy, which is a good model for the faintest galax- 
ies in our sample, then derive approximate scaling relations for 
the dependence of N on M m in nuclei with shallower power-law 
indices. 
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FIG. 6. — Top panels: dependence of consumption rate on M, in NGC 4551 
(power-law galaxy) and NGC 4168 (core galaxy). Dashed lines are the ap- 
proximate Mm -dependences derived in §5.2. Bottom panels: M.-dependence 
of three characteristic radii: r,, the tidal disruption radius; r cru , the radius di- 
viding the diffusive- and full-loss-cone regimes; and r/,, the black hole's radius 
of influence, rj, is the break radius of the luminosity profile. 



5.1. The Singular Isothermal Sphere 

Faint galaxies like M32 have the greatest consumption rates. 
These galaxies also have steep power-law nuclear density pro- 
files, p ~ r" 7 , 7 fa 2. Since r crlt is generally less than rt, in our 
galaxies (cf. Figure [6jl, we can approximate the stellar density 
profile as a single power law when computing the flux. Here we 
consider the singular isothermal sphere (SIS), 7 = 2. The SIS 
density profile and potential are 



P(r) = 



2irGr 2 



ip*(r)=-2a 2 In 



GM. 



(16) 



where a is the ID stellar velocity dispersion, independent of ra- 
dius for r>r/ t . The potential due to the stars has been normal- 
ized to zero at r = r%, and ip(r) = ip*(r) + GM m / r. The isotropic 
distribution function describing the stars is 



1 



dp dip 
V87r 2 m i de J dip y/e-ip 



rla 3 \ m* 



(17a) 
(17b) 



*(£*) = 



y/2 r* L 2 (u)[2+L(u)] dip* 
4^ Zoo [\+L{u)f V^W'' 



(17c) 



L(u) is the Lambert function (also called the W function) de- 
fined implicitly via u = Le L . The superscript "*" denotes dimen- 
sionless quantities and the units of mass and velocity are 



[M]=M„ [V] = a 



(18) 



with G = 1. The dimensionless function ^(e*) that characterizes 
the deflection amplitude per orbital period is 



3V2 \M.J ip*(r,)-e* \r h 
and the loss cone flux is 

2567T 4 InA 

where 

h*(e*) = h\(e*) + h* 2 (e*) + h*(e*), 
h*(s*) = 2 



(19) 



(20) 



g*(e*')de*' 



h* 2 (e*) = 3 



h*{e*) = - 



r * (E * ) dr*r* 2 



r(r*)-e*J £ , 



(21a) 
(21b) 



de*'y/W)-e*'g*(e*'), (21c) 



dr*r* 



[iP*(r*)-e*Y 



de*' [ip*(r*)-e*'] 3/2 g*(m'$) 




FIG. 7. — Dimensionless functions g*(e"), h*(e*) that characterize the 
phas e-space density and angular momentum diffusion coefficient (equations 
I17cl21) in a singular isothermal sphere galaxy. 

Note that the functions g*(e*) and h*(e*) are determined 
uniquely in these dimensionless units. These functions are plot- 
ted in Figure [7] 
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FIG. 8. — Dimensionless loss-cone flux T* (equation 1201 and q (equation 
1191 as functions of energy in a singular isothermal sphere galaxy, for various 
values of Mm . The Mm - a relation (equation^ was used to relate a to Mm . 
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FIG. 9. — Consumption rate as a function of M. in singular isothermal 
sphere nuclei, under two assumptions about cr(Mm). Dashed lines: equation 
(29). 



The function Ro(e) that defines the edge of the loss cone is 
given by equation (TT2l . with 



Ric(e) = 2[ - 
Jh 



(2 + r* -1 ), 



-U-£*)/2/4 . 

(22) 

r c (e) is the radius of a circular orbit of energy e. 

If we set A = 0.4M./m* (Spi tzer & Hart 1971| l, the dimen- 
sionless flux T*{e*) is determined by the two parameters 



M. rh 

m* ' r t 



(23) 



Adopting equation for r,, we can write the second of these 
two parameters as 



r h 29 (M % 



V 



2/3 



/;/,„ 



:21.5('^ 



2/3 



2/3 



100 km s" 1 



M P 



(24a) 



— l(24b) 



with = Gm i /2<j 2 r !t the Safronov number; r\ has been set 
to 0.844. Given values for m* and r*, the two parameters that 
specify T(E) are then M. and a. 

Figure [S^ shows T*{e*) for various values of M.; a was 
computed from M, via the M, - a relation Q. The flux ex- 
hibits a mild maximum at s ~ £/, and falls off slowly toward 
large (bound) energies. The function q(e*) is shown in Figure 
ISJ). As M m is reduced, more and more of the nucleus lies within 
the full-loss-cone regime, q 3> 1. 



Figure [9] shows the consumption rate N = J T(E)dE as a 
function of M, under two assumptions about the relation of a to 
Mm'- cr = 100 km s -1 ; or a determined from the M, - a relation 
0. For fixed cr, Figure|9]shows that N ~ M~ l , while allowing a 
to vary with M. implies a weaker (but still inverse) dependence 
of N onM.. 

The scaling of the consumption rate with M. and a can be 
derived in a straightforward way. Figure [8] shows that over a 
wide range of M. values, most of the flux comes from e > e/,. 
In this energy interval, tp*(r*) « GM,/r* and 



V2ir 3 



A/2 



h*(e*) : 



5V2 
24^ ( 



(25) 



h\, i.e. most 



the latter expression makes use of the fact that h 
of the flux comes from scattering by stars with energies greater 
than that of the test star. Thus 



and Ric w A(r, / r/,)e*. The dimensionless flux is 



160 In A 

9\/2tt £ 



A.^lnA 
9 

B = —, 
4r t 



*l/2 
111* 



A + £ * 2 lnf- 



(26) 

(27a) 

(27b) 
(27c) 



Ignoring the weak E-dependence of the logarithmic term and 
taking r/, jr t from equation d24bi . we find 



N = J HE)dE oc A"'/ 4 oc ct 7 / 2 m; 11/12 . 



(28) 
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After some experimentation, we found that the following, 
slightly different scaling: 



/V«7.1 x 10"V 



70 km s" 



7/2 



V10 6 A1 C . 



(29) 



provides a better fit to the exact, numerically-computed feeding 
rates over the relevant range in M, (Figure |9j. The normal- 
ization constant in equation \29\ was chosen to reproduce N 
exactly for a = 70 km s" 1 , M. = 1 6 M. © . The consumption rate 
in equation J29l > was derived assuming stars of solar mass and 

radius and scales as m±^rJ A . 

5.2. Shallower Power-Law Profiles 

The bright galaxies in our sample have nuclear density pro- 
files with shallower power-law indices, p ~ r" 7 , 7 < 2. We 
calculate the dependence of N on M m in these galaxies using a 
more approximate approach. 

The flux of stars into the loss cone, equation 1 11 PL can be 
written 



HE) = 



T max {E) = 4TT 2 P(E)J 2 (E)JI(E)f(E) 
» HE)N(E) 



(30a) 

(30b) 
(30c) 



with N(E) the number of stars per unit energy interval; equa- 
tion ( l30cl assumes that P(E,J) ps P(E). Now \i = 2r 2 (Av 2 ) /J 2 , 
a function both of r and E, and its orbit-averaged value JI can be 
interpreted as the time-averaged inverse of the relaxation time 
Tr for orbits of energy E. Hence 



HE): 



N(E) 



1 



T R (E) InR^(E) 



(31) 



Above some energy E cr j t , Rq 
creasing E, while for E < E a 
varying function of order unity. Hence 



Ri c e q falls off rapidly with in- 

.-1 



R ps R !c and \nR l is a slowly 



N ! 



Ecu 



ME) 
Tr(EY 
N(r < r cril ) 

Tu(/crit) 



dE 



(32a) 



(32b) 



where &(r cr i,) ps E crlt . These expressions correspond physi- 
cally to the fact that the time to scatter into the loss cone is com- 
parable to Tr for all E < E cr it- Frank & Rees (1977, hereafter 
FR) used an equation similar to ( I32bl > to estimate feeding rates 
in nuclei with constant-density cores. We repeat their analysis 
here, for black holes in nuclei with arbitrary density slopes: 



P(r) = po 



(33) 



For Tr we take the Spitzer & Harm (1958) reference time: 

V2a 3 (r) 



T R {r)-. 



TtG 2 m+p(r)\n A 
Since r cn( < r/, (Figure|6|l, we can write a 2 (r) w GM m /r and 



(34) 



N 



■ N(r < r crU ) 



T R (r crit ) 



9/2-27 



(35) 



Following FR, we define r crif to be the radius above which en- 
counters can scatter stars into or out of the loss cone in a single 
orbital period; at this radius, q w 1. At r = r cr ,r, the angular size 
of the loss cone, 8i c , is comparable to the angle Oj by which 
a star is deflected in a single period; taking account of gravi- 
tational focussing, 9{ c 1=3 \/r,/r. We adopt equation ^ for r t 

P/T r 



and write 9 2 d 



r 3 /GM,/T r ; the square root depen- 
dence of 8 C { on P reflects the fact that entry into the loss cone is 
a diffusive process. Setting Q\ c = 64 then gives 



r crit 
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= Constant x (In A)" 1 m~ l p^r^M 1 . 12 (36) 



and for fixed (po, ro), 



NocM 6 . 



27-197 
6(4-7) ' 



(37) 



For the SIS, 7 = 2 and we recover <5 = — 11/12 (eauationl28>. 

For 7 < 27/19 = 1.42, equation (|37} gives 5 > and N in- 
creases with increasing M.; for instance, setting 7 = gives the 
constant-density core and oc Mi' 1 . This explains why the tidal 
destruction rates in the "core" galaxies generally increase with 
increasing M % . As M, is increased still further in these galax- 
ies, drops, since r/, > rt and the effective power-law index 
becomes steeper. Figure 6 shows fits of equation (12 8 1 to N(M,) 
for two galaxies. 

6. IMPLICATIONS FOR THE DETECTION OF FLARES 

Black hole mass is observed to correlate tightly 
with bulge velocity dispersion (Ferrarese & Merritt 20001 
Gebh ardt et al. 200"0l . and with bulge mass (Merritt & Ferrarese 200 la> 
and luminosity ( McLure & Dunlop 2002 Erwin, Graham & Caon 2002 1. 
Hence we can convert our scaling relation ( I29> into a net scal- 
ing of N on L. First combining the M, - a relation (0 with 
equation d29l : 



N w 4.2 x 10" 4 yr 



100 km s" 1 



6.5 x lO^yr -1 



/ M. 



V10 6 A1c 



-0.25 



(38a) 
(38b) 



Merritt & Ferrarese (2001a) find that \og lQ M./M hu i ge is 
distributed as a Gaussian with mean -2.91 and dispersion 
0.45; the latter is consistent with being due entirely to mea- 
surement errors in M,. Magorrian et al. (1998) find a 
mean mass-to-light ratio for their galaxy sample of sa 



,0.18 , 



4.9/z (L/lO I() /r 2 £ ) T (h = H /80 km s" 1 Mpc" 1 ). Com- 
bining these relations with equation J38bi gives 



A> fa 2.2 x lO'V"'/! 



-0.25 



10 10 /r 2 £, 



-0.295 



(39) 



MT derived a similar relation (their "toy model", eq. 58) for 
consumption in a power-law nucleus. Correcting for different 
assumed Hubble constants, their relation is 



N : 



> 2.6 x 10" 5 yr _1 /! 2/3 



lO n) h- 2 Cr. 



-0.22 



(40) 
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the different scaling with h results from their use of an effec- 
tive radius-luminosity relation in place of the M, - a relation. 
Our predicted event rate is a factor ~ 7(12) greater than theirs 
at L = 10 10 (10 8 )£q; these differences result primarily from the 
larger value (0.006 vs. 0.001) assumed by MT for (M./M bulge ), 
and secondarily from our steeper dependence of N on L. 

MT derived a total flaring rate for early type galaxies 
and bulges by combining equation (I40i with the Ferguson & 
Sandage (1991) E + SO luminosity function and assuming an 
equal contribution from black holes in bulges. They found a 
rate per unit volume of 6.6 x 10~ 7 yr" 1 Mpc" 3 (Ho = 80). Com- 
paring equations (I40i and ( I39> . we conclude that the downward 
revision in black hole masses implies roughly an order of mag- 
nitude increase in the total event rate, to ~ 10~ 5 yr" 1 Mpc -3 . 

The faintest systems in which there is solid kinematical ev- 
idence for nuclear black holes are M32 and the bulge of the 
Milky Way (L w 1O 9 £ ,M. w lO 65 ^©). However there 
is compelling circumstantial evidence for supermassive black 
holes in fainter systems (e.g. Filippenko & Ho 2003) and 
less compelling evidence for intermediate mass black holes 
(IMBHs) in starburst galaxies and star clusters (van der Marel 
2003 and references therein). Here we consider the conse- 
quences for the overall tidal flaring rate if nuclear black holes 
exist in galaxies fainter than M32. The galaxies in ques- 
tion are the dwarf ellipticals, spheroidal systems fainter than 
My sa -19 (Ferguson & Binggeli 1994 1. In spite of their dis- 
tinct name, dE galaxies have properties that are a smooth con- 
tinuation to lower luminosities of the properties of bright ellip- 
ticals dJerjen & Binggeli 1997]lGraham & Guzman 20031 . The 
dEs are the most numerous type of galaxy in the universe; in 
rich galaxy clusters their numbers appear to diverge at low lu- 
minosities as N(L) ~ L~ l (Ferg uson & Sandage 1991^ . If dE 
galaxies contain nuclear black holes, they would dominate the 
total tidal flaring rate due both to their numbers and to their high 
individual event rates. 

Rather than assume that every dE contains a nuclear black 
hole, we make the more conservative assumption that only the 
nucleated dEs (dEn's) contain black holes. Most of the dEn's 
are too distant for their central luminosity profiles to be re- 
solved (e.g. Stiavelli et al. 2001); one exception from the 
Local Group is NGC 205, in which the deprojected density is 
observed to increase as ~ r~ 2 inward of ~ 1 pc (L. Ferrarese, 
private communication). This is similar to what is seen in the 
other, nucleated spheroidal systems in the Local Group with 
comparable luminosities, namely M32 and the bulges of M3 1 
and M33 (fLauer et al.~i9 98) and the bulge of the Milky Way 
dGenzel et al. 2 003 ). We assume that all dEn's have nuclei with 
a similar structure and that dEn's contain nuclear black holes 
with the same ratio of black hole mass to total luminosity that is 
characteristic of brighter galaxies. We can then apply our scal- 
ing relations, equations d38al bl3*9l. to dEn galaxies. We note in 
passing that the luminosity profiles of the dEns - a steep nu- 
cleus superposed on a shallower background profile - is just 
what is predicted by "adiabatic growth" models for black holes 
dPeebles 19721 1 Young 19 80 1, although there are other possible 
explanations for the origin of the nuclei (e.g. Freeman 1993). 

Van den Bergh (1986) plots the fraction of dEs that are nu- 
cleated in a sample of galaxy clusters observed by Binggeli, 
Sandage & Tammann (1985). He finds a roughly linear relation 
between the nucleated fraction F„ and absolute magnitude: 

F n K -0.2 (M v + 13), -18 < My < -13. (41) 

The nucleated fraction is unity in dEs brighter than My w -18 



and negligible in galaxies fainter than My w — 13. Trentham & 
Tully (2002) derive luminosity functions for the dwarf galaxy 
populations at the centers of six galaxy clusters including the 
Virgo cluster. They fit their data to a Schechter function, 

/ L \ ad+1 

N(M R )dM R =N d (—J e" L/i " dM R (42) 

with Mr the 7?-band absolute magnitude. The normalization 
factor Nd has units of Mpc" 2 and gives the surface density of 
dE galaxies at a distance of 200 kpc from the cluster center. 
Trentham & Tully find faint-end slopes of — 1.5 < ay < — 1, con- 
sistent with earlier determinations (e.g. Sandage, Binggeli & 
Tamman 1985). 

Table 2 gives the tidal flaring rate implied by equations ( I39> . 
J41i and J42I for the centers of the six clusters analyzed by 
Trentham & Tully (2002). We give also an event rate for the 
center of the Coma cluster based on the Seeker & Harris (1996) 
dE luminosity function. The highest events rates, ~ 0.1 yr -1 
Mpc" 2 , are predicted for the centers of the Coma and Virgo 
clusters. 

These predicted event rates could be substantially in- 
creased by including the contribution from the bulges of 
late-type spirals, assuming the latter also contain IMBHs. 
Bulge luminosity profiles are similar to those of dE galaxies 
(Mollenhoff & Heidt 200 fl IBalcells et al. 2003> and often ex- 
hibit distinct nuclei dCarollo et al. 20021 . Balcells, Dominguez- 
Palmero & Graham (2001) present resolved nuclear density 
profiles in a sample of spiral bulges observed with HST; the 
nuclei are well fit by power laws with 1 .5 < 7 < 2.5, similar to 
what is seen in the nuclei of the Local Group dwarves. 

The event rate due to dEn galaxies in the Virgo cluster as a 
whole can be computed using the determination by Ferguson & 
Sandage (1989) of the spatial distribution of the dEn's. They 
find a surface density £(/?) w E e" R/s °, ^0 ~ 0.48 Mpc, more 
centrally concentrated than the distribution of non-nucleated 
dwarves. Using the central density normalization of Trentham 
& Tully (2002), we find a total rate of tidal flaring due to dwarf 
galaxies in Virgo of ~ 0.16 yr" 1 . Assuming Poission statistics, 
the probability of detecting at least one event would be 0.15, 
0.55 and 0.80 after 1, 5 and 10 yr respectively in the Virgo 
cluster alone. While the spatial distribution of the dEn galaxies 
in the Coma cluster has apparently not been determined, we ex- 
pect higher overall rates in Coma than in Virgo due to its greater 
richness. 

Some tidal flaring models f Gurzady an & Ozernoy 1980| 
ICannizzo, Lee, & Goodman 1990) predict that single flares 
should persist for as long as several months or years, and in- 
spection of the light curves of the handful of candidate X-ray 
events ( Komossa & Dahlem 2002 1 suggests decay times of this 
order. Such long decay times would imply a non-trivial prob- 
ability of observing an ongoing disruption event somewhere in 
the Virgo or Coma cluster at any given time. Non-detection of 
X-ray flares in these clusters would constitute robust evidence 
that dE galaxies do not harbor IMBHs. 
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7. CONCLUSIONS 

1 . In most galaxies, the predicted rate of stellar tidal disrup- 
tions varies inversely with assumed black hole mass. This is 
particularly true for galaxies with steep central density profiles, 
which dominate the overall event rate. 

2. The downward revision in black hole masses that followed 
the discovery of the M, - a relation implies a total flaring rate 
per unit volume that is about an order of magnitude higher than 
in earlier studies. 

3. An accurate analytic expression (equation!29> can be de- 
rived that gives the tidal flaring rate as a function of black hole 



mass and stellar velocity dispersion in galaxies with p oc r 
nuclei. 

4. If black holes are present in nucleated spheroids fainter 
than My as — 19, the tidal disruption rate due to dwarf galaxies 
in the Virgo cluster would be of order 0.2 yr -1 . Non-detection 
of flares after a few years of monitoring would argue against the 
existence of intermediate mass black holes in dwarf galaxies. 

We thank B. Binggeli, H. Cohn, A. Graham, M. Milosavlje- 
vic, C. Pryor, and M. Stiavelli for useful discussions. This work 
was supported by NSF grants AST 00-71099 and AST02-0631 
and by NASA grants NAG5-6037 and NAG5-9046. 
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Table 2 

Predicted Event Rates Due to Dwarf Galaxies 



Name 


Rate 




(yr -1 Mpc" 2 ) 


Coma cluster 


0.093 


Virgo cluster 


0.11 


NGC 1407 group 


0.057 


Coma I 


0.032 


Leo group 


0.0063 


NGC 1023 group 


0.027 


Ursa Major group 


0.0079 
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Table 1 
Galaxy Sample 3 



Name 


Profile* 3 


Distance 


lno . ( fi ) 
lu 610v r>J 




a 


g 


r 


Ty 


lno.„f/w / J 
lu 510v * / ' 


lu ol0V • / tyl Q ' 


loc,„ /V d 

lu el0 JV 


lu el0v w » / 1 


iog,„ yv f 






^ivipc ) 


(pc) 










(Mq/Lq) 






(yr ) 




(yr ) 


NGC221 


\ 


0.8 


-0.26 


11.77 


0.98 


1.36 


0.01 


2.27 


8.57 


6.38 


-3.79 


6.32 


-3.78 


NGC224 


\ 


0.8 


0.11 


13.44 


4.72 


0.81 


0.12 


26.1 


9.86 


7.79 


-3.70 


6.13 


-3.56 


NGC596 


\ 


21.2 


2.56 


18.03 


0.76 


1.97 


0.55 


4.16 


10.29 


8.69 


-5.01 


7.69 


-4.52 


NGC1023 


\ 


10.2 


1 .96 


16.17 


4.72 


1.18 


0.78 


5.99 


9.99 


8.55 


-4.46 


8.17 


-4.19 


NGC1172 


\ 


29.8 


2.55 


18.61 


1.52 


1.64 


1.01 


2.57 


10.23 


8.42 


-4.75 


6.90 


-3.24 


NGC1426 


\ 


21.5 


2.23 


17.53 


3.62 


1.35 


0.85 


4.91 


10.07 


8.54 


-4.87 


7.50 


-4.08 


NGC3115 


\ 


8.4 


2.07 


16.17 


1.47 


1.43 


0.78 


7.14 


10.23 


8.61 


-4.19 


8.74 


-4.28 


NGC3377 


\ 
\ 


9.9 


0.64 


12.85 


1.92 


1.33 


0.29 


2.88 


9.81 


7.79 


-4.16 


7.51 


-4.04 


NGC3599 


\ 


20.3 


2.12 


17.58 


13.0 


1.66 


0.79 


2.09 


9.82 


7.91 


-5.24 


6.22 


-4.15 


NGC3605 




20.3 


1.94 


17.25 


9.14 


1.26 


0.67 


4.05 


9.59 


8.10 


-5.17 


6.76 


-4.50 


NGC4239 




15.3 


1.98 


18.37 


14.5 


0.96 


0.65 


3.37 


9.19 


7.49 


-5.57 


5.69 


-4.84 


NGC4387 




15.3 


2.52 


18.89 


3.36 


1.59 


0.72 


5.34 


9.48 


7.99 


-5.13 


6.83 


-4.46 


NGC4434 




15.3 


2.25 
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1.78 


0.70 
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9.52 


7.97 
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-4.16 


NGC4458 




15.3 


0.95 
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0.49 


4.00 
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7.90 


-4.58 


6.78 
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NGC4464 


\ 


15.3 


1 .95 


17.35 


1.64 


1.68 


0.88 


4.82 


9.22 


7.69 
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7.26 


-3.88 


NGC4467 
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15.3 


2.38 
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6.27 
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7.32 
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-3.30 


NGC4478 
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36 4 


2 65 


18 33 


95 


1 50 


14 


7 54 


10 64 
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NGC4365 


n 


22.0 


2.25 


16.77 


2.06 


1.27 


0.15 


8.40 


10.76 


9.46 


-5.29 


8.57 


-5.14 


NGC4472 


n 


15.3 


2.25 


16.66 


2.08 


1.17 


0.04 


9.20 


10.96 


9.42 


-5.15 


8.79 


-5.05 


NGC4486 


n 


15.3 


2.75 


17.86 


2.82 


1.39 


0.25 


17.70 


10.88 


9.56 


-5.35 


9.16 


-5.28 


NGC4486b 


n 


15.3 


1.13 


14.92 


2.78 


1.33 


0.14 


9.85 


8.96 


8.96 


-4.84 


8.17 


-4.43 


NGC4636 


n 


15.3 


2.38 


17.72 


1.64 


1.33 


0.13 


10.40 


10.60 


8.36 


-5.35 


8.07 


-5.37 


NGC4649 


n 


15.3 


2.42 


17.17 


2.00 


1.30 


0.15 


16.20 


10.79 


9.59 


-5.19 


9.19 


-5.12 


NGC4874 


n 


93.3 


3.08 


19.18 


2.33 


1.37 


0.13 


15.00 


11.35 


10.32 


-6.02 


8.77 


-5.91 


NGC4889 


n 


93.3 


2.88 


18.01 


2.61 


1.35 


0.05 


11.20 


11.28 


10.43 


-5.81 


9.20 


-5.69 


NGC5813 


n 


28.3 


2.04 


16.42 


2.15 


1.33 


0.08 


7.10 


10.66 


9.29 


-5.24 


8.27 


-5.10 


NGC6166 


n 


112.5 


3.08 


19.35 


3.32 


0.99 


0.08 


15.60 


11.32 


10.47 


-6.16 


8.84 


-6.14 


NGC524 


n 


23.1 


1.55 


16.02 


1.29 


1.00 


0.00 


14.30 


10.54 


9.47 




8.62 




NGC1316 


n 


17.9 


1 .55 


14.43 


1.16 


1.00 


0.00 


2.56 


11.06 


9.25 




8.36 




NGC1400 


n 


21.5 


1.54 


15.41 


1 .39 


1.32 


0.00 


10.70 


10.36 


9.16 




8.62 




NGC1700 


\ 


35.5 


1.19 


13.95 


0.90 


1.30 


0.00 


4.00 


10.59 


8.97 




8.38 




NGC2636 


\ 


33.5 


1.17 


15.68 


1.84 


1.14 


0.04 


2.97 


9.47 


7.72 




6.52 




NGC2832 


n 


90.2 


2.60 


17.45 


1.84 


1.40 


0.02 


10.90 


11.11 


10.06 




9.05 




NGC2841 


\ 


13.2 


0.92 


14.55 


0.93 


1.02 


0.01 


8.98 


9.88 


8.62 




8.23 




NGC3608 


n 


20.3 


1.44 


15.45 


1.05 


1.33 


0.00 


7.04 


10.27 


8.39 




8.01 




NGC4552 


n 


15.3 


1.68 


15.41 


1.48 


1.30 


0.00 


7.66 


10.35 


8.67 




8.62 




NGC7768 


n 


103.1 


2.30 


16.99 


1.92 


1.21 


0.00 


9.51 


11.10 


9.93 




8.82 





a All parameters except for black hole mass and consumption rate are taken from Faber et al. (1997) and assume Hq = 80 km s 
Mpc" 1 . 

b Profile class: n=core galaxy; \=power law galaxy 
c Black hole mass from Magorrian et al. (1998). 

d Consumption rate based on the Magorrian et al. (1998) black hole mass. 
e Black hole mass from the M, - a relation, equation @. 
Consumption rate based on the M. - a black hole mass. 



